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[l]  We  show  high-resolution  satellite  observations  of 
mountain  wave  events  in  the  stratosphere  above  South 
Georgia  Island  in  the  remote  southern  Atlantic  Ocean  and 
compute  the  wave  momentum  fluxes  for  these  events.  The 
fluxes  are  large,  and  they  imply  important  drag  forces  on  the 
circulation.  Small  island  orography  is  generally  neglected  in 
mountain  wave  parameterizations  used  in  global  climate 
models  because  limited  model  resolution  treats  the  grid  cell 
containing  the  island  as  ocean  rather  than  land.  Our  results 
show  that  satellite  observations  can  be  used  to  quantitatively 
constrain  mountain  wave  momentum  fluxes,  and  they 
suggest  that  mountain  waves  from  island  topography  may 
be  an  important  missing  source  of  drag  on  the  atmospheric 
circulation.  Citation:  Alexander,  M.  J.,  S.  D.  Eckermann, 
D.  Broutman,  and  J.  Ma  (2009),  Momentum  flux  estimates  for 
South  Georgia  Island  mountain  waves  in  the  stratosphere  observed 
via  satellite,  Geophys.  Res.  Lett.,  36,  LI 28 16,  doi:  10.1 029/ 
2009GL038587. 

1.  Introduction 

[2]  Accurate  representation  of  the  stratospheric  circula¬ 
tion  can  have  important  effects  on  seasonal  to  decadal 
climate  forecasts  [ Shaw  and  Shepherd ,  2008;  Deckert  and 
Dameris ,  2008],  and  mountain  wave  breaking  has  a  strong 
effect  on  the  circulation  at  these  levels  [Palmer  et  al. ,  1986; 
McFarlane ,  1987;  Eckermann  and  Preusse ,  1999].  Finite 
computing  resources  limit  spatial  resolution  in  global 
models,  so  drag  forces  from  small-scale  breaking  mountain 
waves  are  parameterized.  The  parameterizations  have  long 
been  recognized  as  critical  for  weather  prediction  [Palmer  et 
al.,  1986;  Webster  et  al.,  2003],  and  now  new  sensitivities  to 
the  drag  are  reported  in  climate  simulations  that  incorporate 
trends  in  greenhouse  gases  and  ozone-depleting  chemicals 
[Sigmond  et  al.,  2008;  Li  et  al.,  2008;  McLandress  and 
Shepherd,  2009].  Mountain  wave  drag  parameterizations 
vary  widely  among  models,  and  global  observational  con¬ 
straints  on  mountain  wave  momentum  fluxes  are  needed  to 
reduce  uncertainties  in  these  important  model  projections. 

2.  AIRS  Observations 

[3]  Global  constraints  can  only  come  from  satellite  mea¬ 
surements.  The  Atmospheric  Infrared  Sounder  (AIRS)  instru¬ 
ment  [Aumann  et  al.,  2003]  on  NASA’s  Aqua  satellite  is 
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unique  in  its  sensitivity  to  wave-like  temperature  anomalies 
at  high  horizontal  resolution,  sufficient  to  image  waves  with 
horizontal  wavelengths  longer  than  40  km  [Alexander  and 
Barnet,  2007].  Satellite  measurements  at  microwave  fre¬ 
quencies  have  a  factor  of  3  coarser  horizontal  sampling 
[Eckermann  et  al.,  2007].  Previous  detection  of  mountain 
waves  from  limb-scanning  satellite  instruments  [Eckermann 
and  Preusse,  1999;  Preusse  et  al.,  2002;  Alexander  et  al., 
2008]  observed  the  wave  vertical  structure  in  fine  detail,  but 
obtained  only  relatively  coarse  sampling  of  the  wave  struc¬ 
ture  in  the  horizontal.  Uncertainties  in  momentum  flux  from 
these  previous  analyses  remain  high  [Em  et  al.,  2004; 
Alexander  et  al.,  2008]. 

[4]  Wing-shaped  wave  patterns  downstream  of  South 
Georgia  Island  (37°W,  54.5°S)  observed  by  AIRS  are  shown 
in  Figure  1.  South  Georgia  is  a  small  island  about  170  km 
long  and  30  km  wide  located  in  the  remote  south  Atlantic, 
roughly  1000  km  to  the  east-southeast  of  the  Falkland 
Islands.  The  highest  point  on  the  island  is  Mount  Paget  at 
2915  m.  The  colors  in  Figure  1  represent  AIRS  brightness 
temperature  anomalies  (J7),  which  are  derived  from  radiances 
in  channel  75  and  that  show  anomalies  in  emission  from  C02 
near  40-km  altitude  in  the  stratosphere  [Alexander  and 
Barnet,  2007]. 

[5]  The  ship  wave-like  patterns  in  Figure  1  persisted  for 
four  days  in  this  event,  observed  6-9  September  2003.  The 
patterns  somewhat  resemble  structure  seen  in  low-level 
clouds  downstream  of  small  volcanic  islands  [Gjevik  and 
Marthinsen,  1978;  Gasso,  2008]  due  to  mountain  waves 
trapped  within  the  lowest  few  kilometers  of  the  atmosphere 
[Simard  and  Peltier,  1982;  Eckermann  et  al.,  2006].  Are  the 
patterns  we  observe  at  40  km  in  the  stratosphere  similarly 
caused  by  trapped  waves  in  a  much  deeper  troposphere- 
stratosphere  duct,  or  are  they  caused  by  vertically  propa¬ 
gating  waves?  If  they  are  vertically  propagating,  they  will 
carry  momentum  fluxes  and  may  subsequently  drive 
changes  in  the  circulation. 

3.  Model  Comparison 

[6]  To  address  this  question,  the  mountain  wave  field  over 
South  Georgia  Island  on  6-9  SEP  2003  was  modeled  from 
0-50  km  altitude  using  a  time-dependent  Fourier-ray  algo¬ 
rithm  [Broutman  et  al.,  2006;  Eckermann  et  al.,  2006; 
Broutman  et  al.,  2009].  Island  topography  was  specified 
using  30"  resolution  United  States  Geological  Survey  eleva¬ 
tion  data  interpolated  to  a  2  x  2  km2  regional  grid.  Wind 
and  temperature  profiles  over  the  island  were  specified 
from  0-50  km  for  each  observation  in  Figure  1  using  the 
closest  six-hourly  NASA  Global  Modeling  and  Assimilation 
Office  (GMAO)  GEOS-4  analyses.  Resulting  wave  solutions 
were  computed  at  1  km  vertical  intervals. 
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Figure  1.  Brightness  temperature  anomalies  /(K)  derived  from  radiances  measured  in  the  AIRS  667.78  cm  1  channel 
that  peaks  near  40-km  altitude.  Areas  where  anomalies  are  smaller  than  3  times  the  measurement  noise  [Alexander  and 
Barnet ,  2007]  appear  white.  Data  in  2003:  6  SEP  0310  UT,  7  SEP  1720  UT,  8  SEP  0300  UT,  9  SEP  0340  UT.  (The  color 
scale  is  saturated  at  ±4.5  K  for  presentation  here,  but  the  maximum  and  minimum  values  range  from  —8  to  ±7  K.) 


[7]  Throughout  6-9  SEP  2003,  winds  over  South  Georgia 
were  ^20  ms-1  at  the  ground,  typically  directed  northeast¬ 
ward  to  eastward.  Strong  persistent  eastward  shear  yielded 
eastward  tropospheric  jets  of  40-60  m  s-1  and  stratopause 
jets  in  excess  of  100  m  s-1  near  40-km  altitude.  In  this 
environment  the  model  continuously  generated  mountain 
waves  that  propagated  rapidly  into  the  upper  stratosphere: 
An  example  for  the  6  SEP  2003  case  is  plotted  in  Figure  2. 
A  surface  Froude  number  close  to  unity  suggests  the 


majority  of  the  pressure  drag  across  the  island  is  expected 
to  result  in  a  mountain  wave  response,  as  implicitly  assumed 
in  the  Fourier-ray  method.  The  wave-induced  temperature 
oscillations  well  reproduce  the  observed  wing-shaped  pattern 
in  Figure  1.  While  the  model  is  capable  of  simulating  wave 
trapping,  the  waves  in  this  wing-shaped  structure  are  freely 
propagating  in  the  vertical  from  the  ground  to  40  km.  This  is 
in  sharp  contrast  to  wing-shaped  tropospheric  cloud  patterns 
observed  downstream  of  islands,  which  are  produced  by 


Figure  2.  Mountain  waves  generated  by  the  Fourier-ray  model.  Only  horizontal  wavelengths  >40  km  are  shown  to  match 
the  resolution  of  the  observations,  (left)  Horizontal  cross-section  of  sensible  temperature  perturbations  Ta  (K)  at  40  km  with 
the  island  at  x  =  y  =  0.  (middle)  Vertical  cross-section  of  Ta  (K)  at  y  =  —80  km.  (right)  Magnitude  of  the  momentum  flux 
vector  (mPa)  computed  from  the  covariance  of  horizontal  (u,  v)  and  vertical  (w)  wind  perturbations  p\((u'w'),  (vW))|  at 
40-km  altitude. 
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Figure  3.  Wave  parameters  computed  for  the  6  SEP  2003  overpass  plotted  in  Cartesian  coordinates  relative  to  the  orbit 
track,  along  an  angle  154°  from  north.  North  is  down  and  to  the  right,  and  S.  Georgia  Island  lies  at  (0,0).  Areas  where 
amplitudes  are  smaller  than  3  times  the  measurement  noise  appear  white.  From  left  to  right:  T\  XH,  Az,  and  the  magnitude  of 
the  momentum  flux. 


trapped  waves.  This  finding  permits  the  use  of  linear  theory 
for  vertically  propagating  waves  in  the  following  analysis  of 
the  AIRS  observations.  The  linear  theory  will  allow  us  to 
convert  satellite  brightness  temperature  amplitudes  to  true 
temperature  amplitudes  for  a  quantitative  comparison  to  the 
model  result,  and  for  computation  of  wave  momentum  flux 
directly  from  the  observations. 


4.  Momentum  Fluxes 

[8]  To  compute  momentum  flux  from  the  AIRS  bright¬ 
ness  temperature  anomalies,  the  full  three  dimensional 
structure  of  the  wave  must  be  determined  [Ern  et  al, 
2004].  Wavelet  analysis  allows  calculation  of  the  spatially 
resolved  brightness  temperature  amplitude  (T7),  horizontal 
wavelength  (A #),  and  propagation  direction  of  waves  in  the 
AIRS  data  [Alexander  and  Barnet ,  2007].  We  derive  these 
properties  using  the  S-transform  [Stockwell  et  al .,  1996], 
and  results  for  the  event  observed  on  6  SEP  0300  UT  are 
shown  in  Figure  3.  Propagation  directions,  visualized  as 
vectors  perpendicular  to  the  lines  of  constant  phase  in 
Figure  1  can  only  be  computed  with  180°  ambiguity,  how¬ 
ever  the  ambiguity  is  easily  broken  with  knowledge  of  the 
background  wind  vector  and  the  assumption  that  mountain 
waves  propagate  in  the  upstream  direction  against  the  wind. 
The  flux  vector  directions  therefore  vary  from  northwestward 
north  of  the  island  to  southwestward  south  of  the  island. 

[9]  Vertical  wavelength  is  also  required  for  the  calcula¬ 
tion  of  momentum  flux.  To  estimate  this  from  the  AIRS 
horizontal  cross-sections,  we  need  linear  theory  and  sup¬ 
plementary  wind  measurements,  which  we  obtain  from  the 
GEOS-4  assimilation.  The  linear  dispersion  relation, 

fN1  \ 

Xz=27r[jj2-kHj  >  (F 

relates  the  vertical  wavelength  to  buoyancy  frequency  (N), 
horizontal  wavenumber  (kH  =  27t/A#),  and  the  component  of 
horizontal  wind  in  the  direction  of  wave  propagation  (U). 
The  wave  propagation  direction  and  horizontal  wavenumber 
are  determined  via  the  wavelet  analysis  of  the  AIRS  image 
[Alexander  and  Barnet ,  2007],  and  we  derive  U  and  N  from 
GEOS-4  analyzed  winds  and  temperatures.  The  vertical 
wavelengths  computed  via  equation  (1)  are  also  shown  in 
Figure  3. 


[10]  Brightness  temperature  anomalies  T\  shown  in 
Figures  1  and  3,  were  computed  from  AIRS  radiance  anoma¬ 
lies  approximated  as  small  perturbations  via  the  Planck 
function  [Alexander  and  Barnet ,  2007].  The  brightness 
temperature  amplitudes  (T7)  are  attenuated  relative  to  true 
atmospheric  temperature  amplitudes  because  of  the  deep 
weighting  functions  associated  with  the  AIRS  channel. 
Attenuation  is  strongest  for  short  vertical  wavelength  waves 
and  weaker  for  very  long  vertical  wavelength  waves.  The 
attenuation  factor  A( Az)  varies  smoothly  from  0  (complete 
attenuation)  at  short  Az  to  1  (no  attenuation)  at  extreme  long 
Az.  For  this  AIRS  channel  and  for  Az  =  20  ±  5  km  typical  of 
these  mountain  waves,  A( Az)  =  0.3  ±0.1  [Hoffmann  and 
Alexander ,  2009].  The  true  temperature  anomalies  Ta  can 
then  be  computed  as  Ta  =  T/A(XZ).  This  gives  true  temper¬ 
ature  amplitudes  Ta  of  1 7  ±  4  K  averaged  over  the  central 
region  of  the  wave  field,  in  close  agreement  with  the  model 
result  in  Figure  2  at  40-km  altitude. 

[11]  With  the  full  three-dimensional  structure  of  the  wave 
event  now  characterized,  we  can  compute  the  vector  momen¬ 
tum  flux  at  the  40-km  altitude  level  via  [Ern  et  al ,  2004], 


(Fx,Fy)=-p 


1  ikf)(g\2 


(2) 


where  ( k ,  l)  is  the  horizontal  wavenumber  vector,  m  =  27t/Az, 
g  =  gravitational  acceleration,  and  f  p  are  temperature  and 
air  density  for  the  large-scale  background  atmosphere 
(excluding  the  waves),  here  defined  by  the  GEOS-4  analysis. 

[12]  The  absolute  value  of  these  fluxes  \{FX,  Ff)\  shown  in 
Figure  3  reaches  very  large  values  ~1000  mPa  in  the 
localized  region  near  the  center  of  the  wave  pattern  where 
the  observed  horizontal  wavelengths  are  shortest.  These 
fluxes  are  again  very  similar  to  the  modelled  fluxes  shown 
in  Figure  2.  For  reference,  zonal-mean  fluxes  for  small-scale 
gravity  waves  are  estimated  to  be  ~  1  - 1 0  mPa  in  the  lower 
stratosphere  [Ern  et  al ,  2004;  Alexander  et  al ,  2008].  Note 
however  that  meaningful  values  of  momentum  flux  should 
be  averaged  over  the  horizontal  wavelength  of  the  waves, 
while  in  Figure  3  we  show  the  fully  resolved  fluxes  at 
^20-km  resolution.  Averaged  over  a  typical  global  climate 
model  grid  scale  of  2.5°  x  2.5°  downstream  of  the  island 
gives  an  average  momentum  flux  of  210  ±  20  mPa,  a  value 
five  times  smaller  than  the  maximum  values,  but  still  large 
compared  to  zonal  means.  The  comparable  average  derived 
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Table  1.  Momentum  Fluxes  Over  South  Georgia  Islanda 


Date-Time 

Max  (mPa) 

2.5°  Avg  (mPa) 

09/06/03-03:10  UT 

1400 

210 

09/07/03-17:20  UT 

640 

80 

09/08/03-03:00  UT 

810 

220 

09/09/03-03:40  UT 

240 

70 

aLocal  maximum  (at  20-km  resolution)  and  2.5°  x  2.5°  averaged. 


from  the  Fourier-ray  model  is  120  mPa.  Table  1  lists  the 
momentum  fluxes  estimated  from  the  four  events  shown  in 
Figure  1,  both  maximum  fluxes,  and  2.5°  x  2.5°  averaged 
fluxes.  Errors  estimated  as  standard  deviations  about  the 
local  means  range  ~  10-20%. 

5.  Discussion  and  Conclusions 

[13]  The  AIRS  observations  suggest  these  mountain 
waves  may  occur  commonly  above  South  Georgia  Island. 
We  examined  additional  AIRS  measurements  from  July  and 
August  when  the  climatological  zonal  winds  are  strong  and 
mountain  waves  can  most  easily  propagate  into  the  strato¬ 
sphere  for  detection  in  AIRS  data.  Mountain  waves 
appeared  in  16%  of  the  AIRS  overpasses  in  July  and  34% 
in  August.  The  AIRS  observations  therefore  suggest  that 
mountain  waves  from  South  Georgia  propagate  to  ^40-km 
altitude  commonly  in  winter. 

[14]  Above  40  km,  the  strong  eastward  winds  in  the 
GEOS-4  analysis  rapidly  weaken  to  ^0-30  m  s-1  at 
60  km  altitude  in  the  local  region  above  the  island  during 
6-9  SEP  2003.  These  high-altitude  shear  reversals  are 
believed  to  be  driven  by  gravity  wave  drag.  Fourier-ray 
mountain  wave  solutions  extending  up  to  60  km  diagnose 
convectively  unstable  waves  above  ^50  km,  which  will 
lead  to  flux  dissipation  and  wave  drag.  Assuming  an  areal- 
and  time-mean  wave  flux  of  80  mPa  (Table  1),  if  just  10% 
of  this  flux  is  uniformly  dissipated  between  50-60  km 
altitude,  these  mountain  waves  exert  a  localized  drag  on 
the  flow  of  ^200  m  s-1  day-1,  more  than  sufficient  to 
drive  the  weakened  circulation  at  these  altitudes.  Compa¬ 
rable  surface  stresses  (^50-100  mPa)  occur  over  southern 
hemisphere  orography  in  climate  model  parameterizations 
[Webster  et  al,  2003]. 

[15]  From  these  observational  estimates  we  conclude  that 
South  Georgia,  and  other  mountainous  small  islands,  are 
recurrent  significant  sources  of  mountain  wave  momentum 
flux  and  climatologically  significant  sources  of  drag  on  the 
circulation  of  the  stratosphere  and  mesosphere.  Existing 
mountain  wave  drag  parameterizations  are  problematic  for 
isolated  islands  like  South  Georgia.  At  low  horizontal  reso¬ 
lutions  typical  of  climate  models,  land-sea  masks  will  gen¬ 
erally  tag  the  large  grid  box  encompassing  South  Georgia  as 
an  ocean  point,  so  the  island’s  mountain  wave  drag  would  be 
neither  resolved  nor  parameterized  by  such  a  model.  In  much 
higher  resolution  models  that  can  identify  the  South  Georgia 
grid  box  as  a  land  point,  the  mean  subgrid-scale  orographic 
variance  used  by  the  parameterization  will  still  be  very  small 
due  to  the  still  large  fraction  of  ocean  around  the  island  at 
the  grid  box  edges,  leading  to  substantial  underestimation 
of  the  mountain  wave  flux  from  the  island.  Although  global 
weather  forecasting  models  are  typically  run  at  much  higher 
horizontal  resolution  than  climate  forecast  models,  isolated 


island  topography  could  remain  a  missing  or  underestimated 
source  of  mountain  wave  momentum  flux  in  global  models 
far  into  the  future. 

[16]  Coarse  resolution  global  models  rely  on  parameter¬ 
ized  mountain  wave  drag  to  correct  biases  in  the  zonal  mean 
wind  and  temperature  structure  at  mid-  to  high-latitudes 
[Palmer  et  al.,  1986;  McFarlane ,  1987;  Austin  et  al , 
2003].  This  has  been  effective  in  the  northern  hemisphere, 
where  there  are  large  continental  landmasses  and  numerous 
mountain  ranges.  However  wind  biases  in  the  models  have 
persisted  at  southern  latitudes.  The  biases  have  improved 
somewhat  with  increasing  resolution,  but  still  remain  today 
l Hamilton ,  2008].  The  present  results  suggest  there  is  a  need 
to  include  mountain  wave  drag  from  isolated  island  topog¬ 
raphy  in  the  parameterizations.  While  the  southern  hemi¬ 
sphere  has  far  fewer  continental  mountain  ranges  than  the 
north,  small  mountainous  islands  are  scattered  across  the 
Southern  Ocean,  many  of  which  appear  to  regularly  launch 
waves  into  the  winter  stratosphere  [Wu  et  al ,  2006].  The 
accumulated  deposition  of  momentum  flux  from  mountain 
waves  generated  by  these  island  sources  may  be  an  important 
missing  source  of  mountain  wave  momentum  flux  and  wave 
drag  in  global  models. 
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